Time Evolution of the Neel State 
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A quasionedimentional spin chain (s = |) is considered as a lattice consisting of two sublattices. 
The attention is paid to the states which are pure spin states of the whole lattice and both sublattices, 
the value of the sublattices' spins being maximum. It is shown that the Neel state can be considered 
as a superposition of such states. The exact equation for this superposition coefficients is developed. 

The possibility of the Neel state to be the eigenstate of a Hamiltonian is discussed. Several model 
Hamiltonians are examined, the well known ones and few novel Hamiltonians being considered. The 
time evolution of the Neel state in different models is studied with the help of Fock-Krylov method. 
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I. INTRODUCTION 

When considering the antiferromagnetic state, one 
usually refers to a model of two interpenetrating equiv- 
alent sublattices. The total lattice is in a state with 
zero z-projection of the total spin, whereas each sublat- 
tice is in a state with maximum spin, the z-projection 
of one sublattice being maximum and the other being 
minimum. This state is known as the Neel state. The 
Neel state attracts with it's simplicity and clearness. 
However, its quantum-mechanical realization and genera- 
tion of the corresponding Hamiltonians encounter serious 
problems. For one, the Neel state wave function is not 
an eigenfunction of the system's total spin—. Hence, the 
non-degenerate eigenstate of any Hamiltonian, which de- 
scribes a system with definite total spin, is not the Neel 
state. Therefore, considering various antiferromagnetic 
models one usually uses Hamiltonians with indefinite to- 
tal spin 2 i &±£&ZJLill2.. At the same time, if a Hamil- 
tonian with definite total spin has degenerate spectrum 
the Neel or Neel-like states could be hidden in the degen- 
erate level functional subspace. However, the structure 
of feasible Hamiltonians with degenerate spectrum is not 
at all evident and it is a problem not only to extract the 
given state from the functional subspace but even to rec- 
ognize whether the state is a member of the functional 
manifold. The approach adopted in the present paper 
is not to consider the Hamiltonian cigenstates and com- 
pare them to the Neel state, but to investigate governed 
by the Hamiltonian time evolution of the Neel state. The 
Hamiltonians where the Neel or Neel-like state is either 
stationary or has a long lifetime can be considered as 
describing the antiferromagnetic systems. 

The quasionedimentional spin chain is considered in 
the paper. First, the pure spin states of the chain are 
discussed, the main attention being paid to the states 
with the maximum spin of each sublattice. It is shown 
in the paper that the Neel state can be expressed as the 
particular linear combination of such states and the ex- 
plicit equations for the linear combination coefficients are 
presented. In addition to the pure Neel state various 



Neel-like states are also considered in this section. 

In subsequent sections the Neel state's time evolution is 
investigated with the help of Fock-Krylov method 11 . Var- 
ious spin-model Hamiltonians are employed here. Among 
them are the well-known ones (Heisenberg model, Ising 
model, XXZ-model) , and several new models with the de- 
generate ground state. The influence of the degeneracy 
on the Neel state's time evolution is examined. Finally 
the separable potential^ was employed to obtain models 
with the Neel state as the ground non-degenerate state. 



II. STATES WITH MAXIMUM SPIN OF EACH 
SUBLATTICE 

The spin chain's states with definite total spin and 
maximum spins of each sublattice are considered. It is 
shown, that particular linear combinations of these states 
possess some of the Neel state's properties and, therefore, 
can be considered as a Neel-like state. One such combi- 
nation is the Neel state itself. 

Consider a quasionedimentional spin chain consisting 
of even number AT = 2n of 1/2 spins k = 1, • • • ,N 
(atomic units are used throughout the paper). The chain 
is divided into two sublattices, designated by the index 
I. The first sublattice (£ = 1) comprises spins with odd 
numbers, the second sublattice (£ — 2) comprises spins 
with even numbers. The spin operators for both sublat- 
tices are 



N/2-i 

£ 

i=0 



(1) 



Because operators S 2 , S 2 i, S 2 2, S z commute with 
each other, they have a common set of orthonormal- 
ized eigenfunctions tp(S, Si, £2, M). Consider particular 
eigenf unctions ips = <p(S, Sq, Sq, 0), where S is in the 
range [0, N/2], M = 0, and sublattices spins have maxi- 
mum value S\ = S2 = N/4 = So. There exists only one 
(apart from the phase factor) function ips for each S. It 
can be seen from the momenta coupling rule 
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^s = £ G(S ,M t , S , -M t \S, O^sImAI-m^ (2) 

where G is a Clebsch-Gordan coefficient, is £-th sub- 
lattice spin z-projection and tp^ M , fSo—Mi are nor " 
malized to unity functions of the first and the second 
sublattices correspondingly (the conventional phase fac- 
tor— is used in the present paper). For any Mi there 

exist only one function ip s ^ M . Therefore, the only pos- 
sibility to obtain another function is to couple Si and 
S2 in reverse order, but, according to the properties of 
the Clebsch-Gordan coefficients, this can change only the 
phase of the function <p$. Thus, there are only N/2 + 1 
functions with a definite total spin and maximum values 
of sublattices' spins. 

In what follows we will consider a normalized to unity 
linear combinations of functions (fs 



produce a Neel-like state, i.e. the state which is not ex- 
actly the Neel state but has properties similar to those of 
the Neel state. These properties are: (i) a total chain's 
magnetization is equal to zero, while sublattices' mag- 
netizations are non-zero, (ii) a mean value of the scalar 
product of spins from the same sublattice is equal to 1/4, 
(iii) a mean value of the scalar product of spins from the 
different sublattices is close to -1/4. 



A. SUBLATTICES MAGNETIZATION 

In this subsection it will be shown that the sublattices 
magnetization in state §3§ may be non-zero. 

Magnetization of a system is usually defined as a 
mean value of a system's total spin z-projection operator. 
Hence, the £-th sublattice magnetization in the state i> 
is 



N/2 

* = c s<Ps, (3) 
s=o 



N/2 N/2 

($\St*\*) = J2 E CsCs'ivslSulvs'). (4) 
s=o S'=0 



It will be shown, that the particular set of the coeffi- 
cients Cs yields the wave function corresponding to the Here Si z is the operator of the £-ih sublattice spin z- 
Neel state. Besides, the set of coefficients can be found to projection. One has 



So 



(y>s\St z \vs>) = G(So,Mi, So, -Mt\S, 0)G(S ,M e , S , -M t \S', 0)(<p% tMt \StzW ,u t ) = 

M e = ~S 



So 



Me [ G(S ,M e , S , -M t \S', 0)G(S ,M e , S , -M e \S, 0) - G(S Q , -M t , S , M e \S', 0)G(S o , -M t , S , M e \S, 0) 

M t =0 



/ A S ° 

(l - (_i)4So-s-s j J2 M e G(S , M t , S 0l -M e \S', 0)G(S , M t , S , -M e \S, 0), 

M e =0 



because according to^, 

G(Sa,M e ,S o ,-M t \S,0) = 
= (-1) 2S °- S G(S , —Mi, So, M e \S, 0). 



B. SCALAR PRODUCTS OF SPINS 

In this subsection it will be shown that in the state 
$ from §5§ the mean value of the scalar product of the 
same sublattice spins has the largest possible value, while 
a mean value of the scalar product of the different sub- 
lattices spins can be close to -1/4. 
First of all, note that 



From the equation (|5|) it follows that the sublattice mag- 
netization will be zero if the state <& contains spin states 
ips of the same parity only. To have non-zero sublattices 
magnetization the state $ must contain spin states ips 
with different parity. The particular state <E>, where all 
coefficients Cg are equal to zero except two C5 and Cs< 

with odd S — S' , is an example of the state with non-zero where is a transposition operator. The mean value of 



= i(2P w - 1), 



(0) 



sublattices magnetization. 



operator ^ is in the range [—3/4; 1/4], because the lower 
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and upper bounds of are —1 and +1. Now consider Sj 
and Sfc belonging to the same sublattice number I, Then 



AT/2 N/2 

($| (%s fe ) |$> = £ £ £7SCr S /<^|(Sj,t*)|^) (7) 

S=0S'=0 

where 



So 

(^|(s„s fe )|^<) = E G(S ,M < ,5o,-M / |5',0)x 

M f = -S 

xG(5 , S , -M,|S, 0)(^ Mf | (%, s fc ) |^ M ^) 

(8) 

due to the orthonormality of wave functions of a comple- 
mentary sublattice 3 — t. 

For the matrix element in the right hand side of ([8]) 
one has 



(^£mJ(%A)I^So,M, 



(9) 



for every Sj,Sfc from the £-th sublattice. It is evident 
when Mi is maximum, i.e. Mi = N/4. However, the 

matrix element mJ > ) I ^So M e ) does not depend 
on Mi as can be seen from the following. One has 



(10) 



W _ 



(11) 

where 



B(S ,Mi) = V(So + M(){S -Mi + l). (12) 

and the square root is positive. At the same time the 
equation 



(sj,s k ),Si_ 



(sjjSfc), (sj- x + Sk x ) — i{s jy + s ky ) 



(13) 



is valid for spins from the same sublattice because the 
scalar product (sj,s k ) commutes with any component of 
the vector-operator Sj + Sfc. Therefore, 



{B(So,Mi)Y ^sIm* \St+St- (sj,s k )\<p ( sl Mt ) = r B ( S Ml yp (^ImJ^ s fc )(g| ~ $1 + &)l¥>So,JW>) = (14) 



(5 (5 + l)-Af, 2 + Af,) / (/) ^ (£) (fl ^ w 



It can be inferred from this equality, that the mean value and one can write 
of the scalar product of spins from the same sublattice 
does not depend on Mi in the state ip s ^ Me • According 



G(s, m, s, —m\S, 0)G(s, m, s, —m\S', 0) = 5ss> , 

m 

I 



(<ps\fa,B k )\<PS>) = E G ( S °> Mf - S °> ®)G(Sq, Mi, S , -M^S', 0X4^1(8,-, s k )Wso,u e 



(t) 

M e ,M' e 

j E G(5 , M e , So, -Mi\S, 0)G(S , M t , S , -M e \S', 0) = 5 SS >\. 



4 

Me 



(15) 
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Hence, for the state $ from ([3]) one has 

($|(s 3 ,s fc )|$) = I]T|C s | 2 = i (16) 

s 

Now, suppose, that a spin j is from the first sublattice 
and a spin k is from the second one. In order to calculate 
the mean value of (s_y, Sfc) in this case, note that from ([6]), 
© and (TT3D it follows, that 

for every Me. Hence, function ip s ^ M is symmetric under 
any transposition of it's variables, and, therefore, func- 
tion tps is symmetric under any transposition of spins 
from the same sublattice. One has 

(tps\(sj,s k )\ips') = 

^ N/2-1 N/2 

= (N/2)2 J2 ^iVsl&p+lMlVS') = (17) 
^ ' ' p=0 q=l 

= (iV/2) 2 <yg|(Si,S 2 )|y>gO- 

Obviously, 

(S 1 ,S 2 ) = ±(S 2 -S*-S*), (18) 
and, therefore, 



state <E>. Contrary to that, the equation (|20p indicates 
that mean values of the scalar products of spins from 
different sublattices essentially depend on coefficients Cs 
values. It is easy to present two sets of coefficients result- 
ing correspondingly in positive and negative mean values. 
The simplest sets are the following, all Cs are equal to 
zero except Cjy/2 = 1 in one set and Co = 1 in the other 
set. However, as it was shown above, in these states the 
sublattices magnetization is equal to zero. Still, there ex- 
ist many states where the mean values of the scalar prod- 
ucts of spins from different sublattices are close to — 1/4 
and the sublattices' magnetization is non-zero. These 
states can be referred to as the Neel-like states. 



C. NEEL STATE 

In this section the coefficients Cs in ([3]) are found to 
make function © exactly coincide with the Neel state 
wave function. 

In the case of a quasionedimentional spin chain there 
are two Neel states with the following spin functions 

^nLi = a{a 1 )l3(a 2 )....a{a N . 1 )(3{a N ), (21) 



^Neel = P{oi)a{a 2 )....l3{cj N - X )a{a N ) (22) 

where <7j is a spin variable of an i-th spin. 

Let us show, that the Neel state can be represented in 
the form of a linear combination © with the following 
coefficients: 



Sss> 



(<Ps\(Pj,8k)\<PS>) = 

5(5 + 1) JV(JV/4 + l) 



2(iV/2) 2 



4(N/2y 



(19) 



= 6, 



ss> 



1 _ 1_ 25(5+1) 
' 4 ~ N + N 2 



We should note, that this result does not depend on to- 
tal spin projection's value and is also true for functions, 
corresponding to the maximum sublattices' spins and a 
non-zero total spin projection. 

Thus, when spins j and k belong to different sublattices 
one can write 



(1) 

Neel y 
5=0 



^ I9S(N) 



•<ps; 



N/2 



<4 2 L = £(-i) 



(N/2-S) 



S=0 



9s(N) 



K 



where gs{N) are Wigner coefficients 



9s(N) 



(25+ l)N\ 



{N/2 + S + l)!(iV/2- sy: 



(23) 



(ps. (24) 



(25) 



($l(Sj,Sfc)|*) 



1 

N 



s 



25(5+1) 
N~ 2 



(20) 

Equation (fl"6|) shows, that mean values of the scalar 
products of spins from the same sublattice are always 
positive and have the largest possible value indepen- 
dently on the particular choice of coefficients Cs in the 



and 



N/2 

s=o 



Nl 



(JV/2)!(JV/2)! 



(26) 



This statement is easy to prove. According to the mo- 
menta coupling rule © one has 
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(25 + 1) 



sy.(2s + s + iy: 



E 



(-l) z (5 - Me)\(S + Mf)!(2S - S)\S\S\ \ ( i) (2) 
z!(25 - 5 - z)![(5 - M £ - z)\\ 2 [(S -S + M t + zW j ' Ps °< M '' P8 °<- 1 



(27) 



where z runs over all integers, for which the factorial 
argument is positive, and (p s ^ Mg is wave function of the 
i— th sublattice. 

Consider a scalar product {w^\ e i\<{> s) ■ Calculating it 
with the help of equation (f27f , one can see that because 
of the relation 



VNeel 



(1) (2) 
^Sq^Sq^ Sq. — Sq 



and functions (p s ^ M orthonormality the only non-zero 
term in the sum corresponds to M e — 5 — N/4. Hence, 
the only possible value of z is z = 0, and one has 



(i>Neel\<Ps) = 



(25+1) 



(N/2 + S + l)\(N/2-S)\ 



9s(N) 



K 

Thus, the equation (|23|) is proved (one can easily see that 
the right hand side of (|2"5|) is normalized to unity). 

The equation can be proved similarly. Note, that 
the only possible value of z in this case is z — N/2 — S. 



III. FOCK-KRYLOV METHOD IMPLICATION 

In this section the evolution of the Neel state in dif- 
ferent model systems is examined. Consider the system 
with Hamiltonian H having the discrete and continuous 
parts of the spectrum 



The initial state will decay if the distribution W(E) has 
only continuous part v{E). In this case the probability 
to find the system at a time t in the initial state 



L(t) = < $(0)|$(t) > 



tends to zero when t tends to infinity. Here $(£) is the 
solution of the Schrodinger equation with the initial state 
$(0). The systems considered in the present paper have 
only discrete energy spectrum and therefore any initial 
state does not decay. Taking a particular state for the 
initial state $(0) one can see that L(t) oscillates (peri- 
odically or aperiodically) between two values L m i n and 
L m ax- If L m i n > then at any time there is the non zero 
probability to find the system in the state $(0). More- 
over, in any case the averaged in time probability to find 
the system in the state $(0) is not equal to zero. How- 
ever, considering a property 



F(t) =< $(t)|F|$(t) > 

one can see that the deviation of F(t) from its initial 
value F(0) is only partially determined by L(t) and, 
therefore, the F(t) behavior must not be similar to that 
of L(t). Hence, the time evolution of the property of 
interest F(t) should be examined separately. 

In the present paper the initial state of a system is 
taken in the form 



$(0) = Pn^n 



(28) 



H(x)ip n (x) = E k ip n (x), 

H{x)ij){E,x) = Eip(E,x). 

According to the Fock-Krylov theorem^, the time evo- 
lution of the system is determined by the initial state's 
$(0) energy distribution W(E) 



where P n are some numerical coefficients. The particular 
case of (|2"5|) is the state ([3]). In calculations only one of 

two Neel states (|2H and |2"2"|) . namely ^l el , was taken 
for the initial state because the results for the other Neel 
state will be the same. However, instead one probability 
L(t), two probabilities were considered: the probability 



w(E) = v n 6(E - E k ) + v(E), 



V„ = \< -0n|$(O) > 



(E) = <^(£)|*(0) > 



Neeli 



(29) 



to find the system in the initial state and the probability 



L 2 (t) = m)\^L) 



(30) 
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to find the system in the other Neel state. Moreover, as 
it was discussed before, we are interested in states which 
have non-zero magnetization of both sublattices. There- 
fore, we also calculated the sublattices magnetizations in 
the state 



S iz (t) = (*(i)|ft»|*(«)>, 



(31) 



where I is the sublattice's number. In the present pa- 
per only the systems with conserved z-projection of the 
total spin are considered and in the initial states the z- 
projection of the total spin is zero. Therefore, only the 
results for one sublattice magnetization are sufficient, the 
magnetization of the other sublattice being of the oppo- 
site sign. 



IV. NUMERICAL RESULTS AND DISCUSSION 



Let us consider, first, the Ising model 



H 



N-l 

J E 



>j,zBj+l,z- 



(32) 



which is anisotropic because its Hamiltonian does not 
commute with S 2 operator. We consider the case J > 0. 
The ground state in this case is doubly degenerate and 
both Neel states (|2"Tj) and ([22| belong to it. In this case 

the Neel state VwL/ ener gy distribution w(E) is the 5- 
function, the Neel state is stationary, and 



not stationary and the probabilities L±, L,2 1 and magne- 
tization S\ z depend on time. Their time dependences are 
shown in figures Q] and [5] for two variants of the Heisen- 
berg model. In figure [1] the Heiscnberg model with only 
the nearest neighbors interaction is shown, whereas the 
figure[5]corresponds to the case with the second neighbors 
interaction included, the latter being somewhat weaker 
than the nearest neighbors interaction. Although the 
Neel state in this model does not exactly decay (the spec- 
trum is pure discrete), one can see from figures [1] and [2] 
that the probability to find the system in the Neel state 
quickly drops down and afterwards its time average re- 
mains small. Besides, the time average of the sublattice 
magnetization is almost zero. 
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Lx(i) = l; L 2 (i) = 0; Si„(t) = ^. 
However, due to the degeneracy, the linear combination 



V2 

also can be considered as the ground state, but in this 
state, contrary to the Neel one, the magnetization of each 
sublattice is equal to zero. 

Second, we consider the Heisenberg model 



N-l 



N-2 

E 1 

3=1 



(33) 



which is isotropic as its Hamiltonian commutes with S 2 
operator. The Neel state is not the eigenfunction of this 
Hamiltonian and the energy distribution w(E) has some 
width. The width depends on the particular values of 
Ji , J 2 parameters and in case | J\ \ > \ J2 1 > it could be 
rather large, its maximum being shifted off the ground 
state energy. The Neel state in the Heisenberg model is 



FIG. 1: Time evolution of the Neel state in the system de- 
scribed by Heisenberg Hamiltonian H = Ji Ylj=i (%>%+i) + 
j,Sj + 2) with Ji — 1 and J2 = 0. 



Next, we consider the well-known XXZ-model with 
Hamiltonian 



'N—l N-l 

H = A ( E Sj,xSj+l, x + S j,V S j+l, 
. 3=1 3 = 1 

N-l 

+W Sj,zSj+l,z 
3=1 



(34) 



This model has two parameters A and W and the consid- 
ered above Ising and Heisenberg models Hamiltonians are 
particular cases of (|34[) with A = and A = W respec- 
tively. Therefore, to consider the XXZ-model we selected 
the particular anisotropic case A = 1, W — 5 which is far 
from both Ising and Heisenberg cases discussed above. In 
the selected case the Neel state energy distribution has 
two main weights V\ and V2, which correspond to the 
ground and first excited states. Weights v\ and V2 are 
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FIG. 2: Time evolution of the Neel state in the system de- 
scribed by Heisenberg Hamiltonian H = Ji Ylj=i (%j + 
■h Ef = "i 2 (s J ,s J+2 ) with Ji = 1, J a = 0.25. 



approximately equal to each other and they are at least 
one order of magnitude greater than the weight of any 
other state. Such energy distribution results in particu- 
lar time dependences of L%, L2 and S\ z shown in figure[3] 
The main feature is a harmonic oscillation of L\(t) and 
L 2 (t) superimposed by a small amplitude and high fre- 
quency oscillations. The L rnax of the harmonic oscillation 
is close to 1, and L min is close to 0. The L\(t) would ex- 
perience exactly harmonic oscillation if all weights other 
than v\ and V2 would be exactly zero. In such oscilla- 
tion the frequency is defined by the energy gap between 
the ground and excited states, the L max value is equal 
to 1, and L m i n value is defined by the weights differ- 
ence, being zero when weights are exactly equal. In our 
case weights v± and v 2 are only close to each other and 
therefore, L m i n is small but not equal to zero. The sim- 
ilar oscillation (with opposite phase) of L%(£) indicates 
that the energy distribution of another Neel state 4> N eei 

is similar to that of V'tvLz ^ u t w ith interchanged v% and 
1)2 . From this we can conclude that both ground and first 
excited states are linear combinations of two Neel states 
with approximately equal weights and neither of them is 
close to a single Neel state. This conclusion is confirmed 
by the results of numerical calculations and is in accord 
with the time behavior of the sublattice magnetization 
shown in figure [3J From this figure it is seen that the 
time average of the sublattice magnetization is close to 
zero. The result of numerical calculations is —0.02. Be- 
sides, the sublattice magnetizations in the ground and 
the first excited states in the case considered are exactly 
zero. Hence, although the time average of the probabil- 
ity to find the system in the initial Neel state is essential 



(0.39 according to the numerical calculations), the mag- 
netic properties of the system do not correspond to the 
Neel state. 
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FIG. 3: Time evolution of the Neel state in the XXZ-model. 

The next model was constructed to satisfy two con- 
ditions. First, the system was chosen to be isotropic. 
Hence, the system Hamiltonian must commute with the 
total spin operator. Second, the Neel state energy dis- 
tribution in the system must have the smallest possible 
number of components. For this the system Hamiltonian 
must commute with sublattices square spin operators. 
Indeed, in this case equations (f23|) and (|24f show that 
only states with maximum spins of each sublattices con- 
tribute to the Neel state energy distribution. Based on 
these considerations, the system Hamiltonian was chosen 
in the form 

JV-2 

H = A (%.%+2) + WS 2 + V (Sf + 3|) . (35) 

with three parameters A, V, and W. Considering the 
Neel state we see, that all functions ips from ([2]) onto 
which the Neel state is exactly expanded are eigenfunc- 
tions of this Hamiltonian. At the same time, from the 
equation (fT5|) it follows that the matrix of the first term 
in (|35p , calculated on tp$ , is proportional to the unit ma- 
trix. Besides, all ips are the eigenf unctions of Sf and 
S\ with the same eigenvalue So (So + !)■ Hence, the ma- 
trix of the third term in (f55]l . calculated on ips, is also 
proportional to the unit matrix. Consequently, parame- 
ters A and V define only the position of the Neel state 
distribution w(E) in the operator H spectrum, while 
the parameter W defines the energy distribution width 
AE = WN(N + 2)/4. The Neel state evolution in the 
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considered model is periodic (however not harmonic), be- 
cause energy differences between all tpg energy levels are 
commensurable (AE jk = W(Sj(Sj + 1) - S k (S k + !)))■ 
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FIG. 4: Time evolution of the Neel state with the Haniiltonian 



// = ^Ef=i 2 (%,s j+2 ) + ws* + {v-w) (si + si 



In figurc[4]the time evolution of the Neel state is shown 
for the particular case of Hamiltonian with A = — 5, 
W — 1 V — — 7. The parameters were chosen to place 
the energy levels of all eigenstates ips, which contribute 
to the Neel state, below the rest of the spectrum. There 
are four main weights v k in the energy distribution in 
this case. Therefore, the oscillation shape is more com- 
plicated than in the XXZ model with two main weights, 
and the maximum value of Lzit) is well below 1. The 
time average of the probability to find system in the ini- 
tial Neel state is about 0.29. However, the time average of 
the sublattice magnetization is small (about 0.001). Be- 
sides, the sublattice magnetization in the ground state is 
zero. Therefore, similarly to the XXZ model, the mag- 
netic properties of the system do not correspond to the 
Neel state. 

Next model was made to have the degenerate ground 
state with all and only functions ips belonging to it. In 
this case the Hamiltonian, similarly to (|35[) . must com- 
mute with the total spin operator and with the sublat- 
tices spin square operators. In addition, the ground state 
must correspond to the maximum spins of sublattices and 
to the zero z-projection of the total spin. Besides, all 
other states must have higher energy than the ground 
state. The simplest Hamiltonian satisfying these condi- 
tions is 



H = -{Si + S%) + Si 

Similarly to the Ising model, the Neel state V'veeZ 
stationary state and in this state 



(36) 



is the 



Li(t) = 1; L 2 (t) = 0; 



In this model, contrary to the Ising model, any Neel-like 
state is also the stationary state. However, for any Neel 
or Neel-like state the partner state with the opposite sub- 
lattices magnetization also belongs to the ground state. 
Therefore, their linear combination with zero sublattices 
magnetization is also the ground stationary state. 
Finally, the simple model with Hamiltonian 



H = S2 + Y,\^s) (A-S(S+l))(<p s \ 



(37) 



A < 0. 



was considered where functions (ps are defined by (|2|). 
This model shows, that the time evolution of the Neel 
state helps us to recognize the existence of the stationary 
Neel-like state in the system where the Neel state itself is 
not stationary. The peculiarity of the model is that the 
ground state is degenerate and the corresponding func- 
tional subspace is spanned by some (not all) functions 
(ps. Therefore, the Neel-like states are eigenstates of op- 
erator (|3"7)) and the Neel state itself is not. The time 
evolution of the Neel state in this model is shown in fig- 
ure [5l From this picture one can see that the probability 
Li(t) to find the system in the initial Neel state oscillates 
around some value close to 1 {L\ time average is equal 
to 0.97 according to the numerical calculations) and the 
sublattice magnetization oscillates around a value close 
to the maximum (time average is equal to 1.94), the men- 
tioned values being corresponded to the Neel-like state. 
Hence, the behavior of the Neel state time evolution, sim- 
ilar to that shown in figure indicates that in the sys- 
tem there exists stationary Neel-like state. Note, that 
in this system, similarly to all systems with degenerate 
ground state discussed above, besides the particular Neel- 
like state, there exists the stationary Neel-like state with 
the same energy and op 

In all discussed above models the Neel state cannot 
be considered as stable. Indeed, in all models, except 
Ising model and model with Hamiltonian (|36|) . the Neel 
state is not the eigenstate. In the Ising model and the 
model (|36j) the Neel state is the eigenstate but it belongs 
to the degenerate level with the other Neel state also 
belonging to this level. In this case even if a certain Neel 
state was generated at a time t = it will be destroyed 
by small external perturbations which are always present 
in real circumstances. To be stable the Neel state must 
be the ground non-degenerate state. However, it occur 
Hamiltonian can be transformed to have the particular 
Neel state as the ground non-degenerate eigenstate. To 
show this we will use the results of the papei— , where 
the following Hamiltonian transformation was developed 
to get a desirable set of states for the lowest eigenstates. 
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FIG. 5: Time evolution of the Neel state with the Hamiltonian 



Consider a Hamiltonian H° with eigcnfunctions ipf and 
eigenvalues E° 



V n = Ex-E°, V 22 = \t\ 2 < /2I/3 >, 



V23 — V32 — —t 



t = 



1 



l+</i|/a> 



All other matrix elements Vij are equal to zero, so the 
separable potential (|38|) has only four non-zero compo- 
nents. 

Using these results we can make the Neel state ^ ^eel 
be the ground non-degenerate state. For this we take the 
ground state eigenfunction and eigenvalue of the model 
for ipi and E®. When the model's ground state eigenvalue 
is degenerate we take the closest to ^^iei eigenstate for 

In the Ising model with J > and in the model with 
Hamiltonian ([55]) the Neel state is the ground state eigen- 
function of the Hamiltonian, so tpi = tp® 



f /(!) 
Case /i = V]veei> 



iPnLi- Inthis 



2E\i> 



o„,,(i) 



JVeei' 



t = 1/2, 



and Vu = Ei—E°, V 22 = £?/4, V 23 = V 32 = -1/4. Tak- 
ing the sum of four components one can see that three 
components cancel each other and the potential (f3"8")) in 
this case becomes simply the energy level shift operator 



#V° = Ef^. 

Consider also an arbitrary set of n orthonormal func- 
tions ipk and an arbitrary set of real numbers Ek (fc = 
1, • • • , n). Ir>i^ the potential in the form of the separable 
potential 

3n 

V - E \fi)ViiUi I. (38) 

was found so that the n lowest eigenf unctions and eigen- 
values of the Hamiltonian H' = H° + V 



V=\^li>(Ei-E° 1 )<^li\ 

In all other models the potential has four non- 

zero components. The values of matrix elements Vij 
calculated for all considered above models, except ([52^1 
and (|36p are given in the table fl] 



Model 


Parameters Values 


V u 


V 22 


V23 


([33]) 


Ji = 1; J 2 = 


-1.5 


-3.49 


-0.68 


If33| 


Ji = 1; Ji = 0.25 


-1.0 


-5.40 


-1.74 


l[M| 


4= 1; W = 5 


0.0 


-49.60 


-2.90 


([35]) 


A = -5; W = 1; V = -7 


-3.0 


-108.90 


-0.61 




S fe = 3; A = -10 


-1.0 


-9.93 


-0.50 



H'i/jI = E'tipl TABLE I: Non -zero matrix elements of the separable poten- 

tial |38] for different models. 



are exactly ipi and E^ 



tpi =ipi, El = Ei, i = 1, 2, • • • , n 
while the rest of the H° spectrum is not affected 

&i = E%> i = n + l,---. 

The eigenfunctions ip' t with i > n differ, of course, from 
ip® because former are orthogonal to ipk, while later are 
orthogonal to ip^., k = 1, • • • , n. 

When only one eigenstate is to be transformed (n = 1) 
the functions fi and the matrix Vij, i,j = 1,2,3 can be 
chosen as follows 



CONCLUSIONS 

Considering the quasionedimensional spin chain the 
Neel state was analyzed and its time evolution in var- 
ious models was examined. It was shown that in those 
models, where the Neel state is not the eigenstate of the 
Hamiltonian and, therefore, is not stationary, the time 
averaged magnetic properties of this state do not corre- 
spond to that commonly attributed to the Neel state. In 
those models, where the particular Neel state is an eigen- 
state of the Hamiltonian, and this eigenstate is degener- 
ate, the Neel state is stationary, but another Neel state 
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